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Abstract 

Periodically driven dynamics of a particle moving in the field Coulomb plus confining 
potential is treated for one and three dimensional cases. Critical value of the external field 
strength at which chaotization will occur is evaluated analytically based on the resonance 
overlap criterion. The analysis of the phase-space dynamics is presented. 

1 Introduction 

Periodically driven dynamics has been subject of extensive research for past few decades both 
in classical and quantum contexts [l]-[6]. An interesting feature of the periodically perturbed 
dynamical system is the chaotization of its motion under certain conditions. This implies 
divergence of the phase space trajectories and unlimited energy growth (in the classical case) 
leading to acceleration. Qualitative and quantitative analysis of the model systems such as 
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kicked rotor shows that there is linear time dependence in the average energy of the periodically 
driven system leads to a diffusion in the phase space [6]. Corresponding quantum studies have 
shown that such a diffusion is suppressed in the quantum case [6] . Therefore periodically driven 
systems are convenient testing ground for the study of dynamical chaos in time-dependent case. 
Great variety of periodically driven systems have been treated from the viewpoint of dynamical 
chaos theory from model systems such as kicked rotor [6] , pendula [2] or billiards [4] to realistic 
systems like atoms and molecules [5, 7, 8, 9, 3, 10] and surface-state electrons in liquid helium 
[11] in the time-periodic fields. Comprehensive theoretical and experimental studies of the 
hydrogen atom in a monochromatic field showed that diffusive excitation of the atom leading 
to ionization during relatively long time can occur. This fact makes hydrogen atom in a 
monochromatic field a convenient system for the theoretical and experimental study of chaos 
in the time-dependent dynamical systems [3, 7, 8, 9]. A theoretical analysis, of the behavior 
of a classical hydrogen atom interacting with monochromatic field, based on the resonance 
overlap criterion [9]- [10], shows that for some critical value of the external field strength e cr , 
the electron enters into chaotic regime of motion, marked by unlimited diffusion along orbits, 
leading to ionization. Experimentally, this phenomenon was first observed by Bayfield and 
Koch [12]. Such an ionization was called chaotic [3, 7] or diffusive [9] ionization. During the 
last three decades chaotic ionization of nonrelativistic atom was investigated by many authors 
theoretically [9, 3, 11] as well as experimentally [7, 8, 12]. 

In this paper we study a model which may be considered as a QCD counterpart of the 
periodically driven hydrogen atom. Namely, we address the problem of regular and chaotic 
motion of a particle bounded in the field of Coulomb plus confining potential in the presence 
of a time-periodic external perturbation. Such system may be used for modelling periodically 
driven quark- ant iquark state, so-called quarkonium. 

Using resonance analysis based on the Chirikov criterion for stochasticity we estimate critical 
values of the external field strength at which quarkonium motion enters into chaotic regime. 

Quarkonium in a monochromatic field can be considered as an analog of the hydrogen atom 
in a monochromatic field, in which Coulomb potential is replaced by Coulomb plus confining 
potential. Quarkonia have been the subject of extensive experimental as well as theoretical 



2 



studies for the last two decades [13]- [14]. In the framework of potential model the description 
of quark motion in hadrons is reduced to solving classical or quantum mechanical equations of 
motion with Coulomb plus confining potential [14]. 

Study of the nonlinear dynamics of hadrons is of importance due to the recent advances made 
in creation of hadronic and quark-gluon matters in the collisions experiments of ultrarelativistic 
heavy ions where creation of hadronic or quark-gluon matter is possible [15]. Quarkonia in 
quark-gluon matter can be considered as a system perturbed by time- dependent force, that 
may lead to chaotization of the quarkonium motion. Indeed, the recent studies on quarkonium 
dynamics showed that regular motion can be expected at a small values of color screening 
mass but the chaotic motion is expected at a large one [16]. Periodically driven quarkonium 
can be also realized in the interaction of mesons with laser fields. This paper is organized as 
follows. In next section we will treat a simple model, a one-dimensional quarkonium under 
a periodic perturbation. In section 3 we extend our treatment to the three-dimensional case. 
Some concluding remarks are presented in section 4. 



2 One-dimensional model 



For simplicity we consider first a one-dimensional model described by a potential 

-- + Xx for x > 



V(x) = < 



oo 



for x < 



where Z = |a s , a s being the effective strong coupling constant and A gives strength of 
the confining potential. As is well known, in the case of the hydrogen atom interacting with 
a monochromatic field, one-dimensional model provides an excellent description of the experi- 
mental chaotization thresholds for real three-dimensional hydrogen atom [7, 11, 3]. The same 
success is to be expected in the case of quarkonium. 

The unperturbed Hamiltonian for the above potential is 

~ 2 Z 



H = ^--- + \x. 

2 x 



(1) 
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We will treat the interaction of the system given by Hamiltonian (1) with the periodic 
external potential of the form 

U(x,t) = excosut. (2) 

with e and u being the field strength and frequency, respectively. Thus the total Hamiltonian 
of the periodically driven quarkonium is 

H = H + U(x,t) (3) 

Formally, the Hamiltonian (1) is equivalent to that of the hydrogen atom in constant ho- 
mogenous electric field. Chaotic dynamics of hydrogen atom in constant electric field under the 
influence of time-periodic field was treated earlier [17, 18]. To explore periodically driven dy- 
namics of our system we introduce so-called action-angle variables and rewrite its Hamiltonian 
in terms of these variables. 

The action variable is defined as [5]: 



n = Li j 2(Ho _ v[x))dx = / > - «>(« - c > ix , (4) 
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where the constants a and c are the turning points of particle and are given by 



H + J Hi + AZ\ H - J Hi + AZ\ 
a = 2^ ' C= 2^ 

Since c < for the action we have 
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- J ^2(H - V(x))dx 
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where 
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here E(k) and K(k) are the elliptic integrals [19] and 
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The angle variable is defined as [5] 



dn 



where 



S = J X p( x ,H(n))dx 
Then in terms of n and x variables the angle can be written as 

9(x, n) = 



1 1 E(k)-K(k) dk 1 E(k)-kVT^¥K(k) dk 
B i X+ x [{X+ x } k In + x {k) kVT^W dn~ l 



(7) 



where k — k(n) is defined by the eq. (6). 

Furthermore, we consider the following two cases: a» 1 and a <C 1. For both cases we 
obtain approximate expression for unperturbed Hamiltonian as a function of angle variable. To 
do this we we solve the eq.(5) with respect to H and use asymptotic estimates for the elliptic 
integrals, E(k) and K(k) [19] for a « 1 and a » 1 that correspond to k « 1 and k m 1, 
respectively. 

For a> 1 we have: 



H = Z 2 An 2/3 ■ 



\ln(AB~ 2 / 3 n 2 / 3 ) 



A 2 n 4/3 



(8) 



with 



A = 



3ttA 3 

2y/2 



Corresponding proper frequency is 

2 



n 



1/3 ^ n 5/3l 



(9) 



For second case, a < 1 we get 



H = 0.5Z 2 (9.7Xn 2 -n- 2 ); 



(10) 



The proper frequency for this Hamiltonian is 



coo = Z 2 (n~ 3 + 9.7n\). 



(11) 
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The Hamiltonian of the perturbed system can be written as 

H = H + e ^2 XkCos(k9 — out), 

with 



(12) 



rZir 

xdn) = - / x(n,9)e tke d9 (13) 
Jo 

being Fourier amplitude of the perturbation. For a< 1 we have 

. AE(n) 1 . 2 irkV\ 
x k (n)^ \k ~2~ ' ^ 

For a » 1 we obtain 

x k (n) = (15) 

As is well known [1, 2, 5], dynamics of a periodically driven system may become chaotic 
if the resonance condition is fulfilled and external field strength exceeds some critical value. 
To estimate this critical field strength, e cr for our system we use Chirkov's resonance overlap 
criterion [1, 5, 9], which can be written as: 

Az/ fc + Az/ fc+ i 

Ti 7^ 77T > 2 A (1°) 

0J (k + 1) -Wo(fc) 

with 

Ai* = (^) 
being the width of the fc-th resonance [1, 5] and 

u' = du /dn. 

From the resonance condition we get 

u {k) - u (k + 1) = — 



k k+1 k(k + l) 
Applying this criterion to our system give by (12) we have for a ^> 1 

0.07ZWA k(k + l) 

6cr ~ n 2 (k + 1) 2 + k 2 X 
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and for a< 1: 
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Table 1 presents the values of the critical field for a « 1 uu, dd, ss, cc and bb quarkonia at 
the following values of parameters: u = 10 9 Hz, a s = 0.1 12 [20], A = 0.2GeV 2 , n = 5; 7; 10. For 
light (u, d, s) quarkonia we use formula (17) and formula (18) for bb- and cc- quarkonia. In the 
Table 2 the critical values of the field strength are displayed for the same values of parameters 
but at a » 1. The huge difference between the corresponding values in these two tables can 
be explained by the fact that these two limit cases make our system (to some extent) equivalent 
to hydrogenlike atom(a << 1) and highly excited confined(by linear potential only) system. 
Therefore the difference between these two cases is caused by the big difference between the 
binding energies. 

In Fig.l the critical field strength is plotted as a function of the action n. Again, three cases 
are considered: periodically driven motion in the field Coulomb potential (A = 0), periodically 
driven motion in the field of Coulomb plus linear potential for the cases a < 1 and for a ^> 1. 

The parameters Z and A are chosen as: Z = 0.15; A = 0.4. As is seen from this figure, 
unlike to periodically driven hydrogen-like atom, where field strength is a decreasing function 
of n, for periodically driven quarkonium field strength increases first up to certain value of n, 
then decreasing occurs. The curve for quarkonium at a « 1 is closer to that for hydrogen-like 
atom, while for a » 1 the difference between the critical fields for quarkonium and hydrogen- 
like atom becomes significant. This can be explained by the fact that for a « 1 the energy of 
the unperturbed quarkonium becomes closer to that for hydrogen-like atom. 

It is well known that the phase-space trajectories of the regular motion lie on tori (so-called 
KAM tori). According to Kolmogorov-Arnold-Moser theorem for sufficiently small fields most 
of the trajectories remain regular. If the value of the external perturbation exceeds some value, 
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which is called the critical field strength, KAM tori start to break down and chaotization of 
the motion will occur [5]. 

In Figs 2-4 the phase-space portraits are plotted and compared for periodically driven pure 
Coulomb potential (hydrogenlike atom), quarkonium system at a <C 1 and quarkonium system 
at a ^> 1 for various values of the ratio e cr /e. 

Fig. 2 presents phase-space portraits the case e cr /e = 0.1; the panels a,b,c are the phase 
space portraits at hydrogenlike atom, quarkonium for a <C 1 and quarkonium for a >> 1, re- 
spectively. Figs. 3 and 4 correspond to the cases when e cr /e = 0.5 and e cr /e = 0.9, respectively. 

The following values of the parameters Z, A were chosen for the plots: Z = 0.15; A = 0.4. 

As is seen from these plots periodically driven motion in the field of Coulomb potential is 
more chaotic than that in the field of Coulomb plus linear potential. In other words confining 
potential makes periodically driven dynamics more regular. Among two limit cases, a < 1 and 
a ^> 1, dynamics for a <C 1 is more chaotic than that for a >> 1, which also confirms that 
confining force suppresses chaotization of motion. 



3 Three-dimensional model 

The Hamiltonian for the three-dimensional model is 

2 r r 

where L is the orbital angular momentum and p r is the radial momentum. 
The action can be expressed in terms of elliptic integrals [21]: 



n 



= J Vdr = j 'J2(E + j -Xr)dr = 



(2Z/3 - L 2 /c + Ec/3)K(k) + E(a - c)/3E(k)+ 

+L 2 (c- 1 -b- 1 )l[({3 2 ,k)]g/V\ (19) 

with a and c being the turning points, K, E, T7 are complete elliptic integrals of the first, 
second and third kind, respectively [19], and the constants are given as 

k 2 = (a-b)/(a-c), 
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(3 = ck 2 /b, 



g = 2/Va - c. 

From eq.(19) the unperturbed Hamiltonian as a function of n can be found approximately for 
E/\ 1, (which corresponds to n > 1 or large quarkonium masses): 



Ho = q\n)^ 



The proper frequency is 



LU 



dH ,2A 2 , 



1/3 



~^ 3n. 



n -V3 _ !^ n -4/3 
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Then the Hamiltonian of the three-dimensional quarkonium in a monochromatic field can 
be written as 

H = H + ea cos (cut) x 



where 



{ — -esin(/9 + 2 ^[xfe sin ip cos k<p + Uk cos ip sin k(j)]}, 



Xk = [ e tWokt xdt, y k = ( e tu >° kt ydt, 

UnkTJ ' yk uakT J y ' 



(20) 



uj kT J ' a u kT . 

and ^ and (j> are the Euler angles. Again, using the resonance overlap criterion (16) in which 
the resonance width is defined by 



a i eTk \ 



where 



r k = \fx 2 k + yl, 
we obtain an estimate for the critical field: 
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This estimate for the critical field assumes that n » 1. If the external field strength has 
the value exceeding e cr , breaking of KAM surfaces in the pase space will occur and and the 
quarkonium diffuses in action and the motion becomes chaotic. 

4 Conclusions 

In this work we explored periodically driven dynamics of a particle bounded in the field of 
Coulomb plus linear potential. Using resonance overlap criterion the estimates of for the critical 
field strength at which dynamics enters into chaotic regime of motion are obtained. The analysis 
of phase-space dynamics by plotting phase space portraits shows that the periodically driven 
dynamics of quarkonium is more regular compared to that for hydrogen-like atom. Also, there 
appears a pick in n-dependence of the critical field strength, that can be explained by the role of 
confining potential. The obtained results are applied for the estimation of critical field strength 
needed for chaotization of various quarkonia. The model considered in this work can be used 
for modelling nonlinear dynamics of quarkonia perturbed by the field of a quark-gluon plasma 
or laser radiation. 

5 Acknowledgements 

This work is supported by the NATO reintegration grant (Ref. No. 981787), by INTAS Young 
Scientist fellowship (Ref. No. 06-1000023-6008) and by a grant of the Uzbek Academy of 
Sciences (FA-F2-084). The work of FCK is supported by NSERCC. 

References 

[1] B.V.Chirikov Phys.Rep., 52, 159, (1979) 
[2] D.F. Escande, Phys. Rep. 121, 167, (1985) 

[3] G.Casati, B.V.Chirikov and D.L.Shepelyansky, Phys.Rep., 154, 77, (1987) 



10 



[4] B. Eckhardt, Phys. Rep. 163, 207, (1988) 

[5] R.Z.Sagdeev, D.A.Usikov, G.M.Zaslavsky, Nonlinear Physics: from pendulum to turbu- 
lence and chaos (Academic Publisher, NY 1988) 

[6] G.M.Izrailev, Phys.Rep., 196 299 (1990) 

[7] R.V.Jensen , S.M.Susskind and M.M.Sanders Phys.Rep. 201, 1 (1991) 

[8] P.M.Koch, K.A.H. van Leeuwen Phys.Rep. 255 (1998) 289 

[9] N.B.Delone, V.P.Krainov and D.L.Shepelyansky, Usp. Fix. Nauk. 140, 335, (1983) 

[10] D.U.Matrasulov, Phys. Rev. A 50 700 (1999) 

[11] R.V.Jensen, Phys. Rev. A, 30, 386, (1984) 

[12] J.E.Bayfield and P.M.Koch, Phys. Rev. Lett, 33, 258, (1974) 

[13] ALEPH Collaboration, R.Barate, et al, Phys.Lett. B, 425 215 (1998). 

[14] S.N.Mukherjee et.al. Phys.Rep. 231 (1993) 203 

[15] NA50 collaboration, Phys. Lett. B 477 28 (2000) 

[16] Jian-zhong Gu et al, Phys.Rev.C, 60 035211 (1999). 

[17] G. P. Berman, G. M. Zaslavskii, and A. R. Kolovskii, Sov. Phys. JETP, 61 925 (1985) 

[18] Mark J. Stevens and Bala Sundaram, Phys. Rev. A, 36, 417, (1987) 

[19] M.A. Abramowitz and LA. Stegun, Handbook of mathematical functions, Nat. Bur. Stand. 
Washington D.C.,1964; 

[20] Particle Data Group, R.M.Barnett et al Phys.Rev.D, 54 1 (1996) 

[21] M.Seetharaman, R.Raghavan and S.S.Vasan J.Phys.A 16 455(1983) 



11 



TABLE 1. The values of the critical field strength for heavy quarkonia. 



No 


Quarkonium 


Quark mass (in MeV) 


Critical field (V/fm) 








n = 5 


n = 7 


n= 10 


1 


cc 


300 


1.215 


0.6163 


0.3008 


2 


bb 


1560 


5.761 • 10 2 


2.901 • 10 2 


1.407 • 10 2 



TABLE 2. The values of the critical field strength for light quarkonia. 



No 


Quarkonium 


Quark mass (in MeV) 


Critical field (V/fm) 


n = 5 


n = 7 


n = 10 


1 


uu 


1 


1.018 • 10 19 


5.192 • 10 18 


2.544 • 10 18 


2 


dd 


2 


8.141 • 10 19 


4.153 • 10 19 


2.035 • 10 19 


3 


ss 


30 


9.158 • 10 22 


4.673 • 10 22 


2.29 • 10 22 
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Figure 1: Critical field strength (in the system of units where quarkonium mass is equal to 1) 
as a function of n, for periodically driven hydrogenlike atom, quarkonium in a ^> 1 case and 
quarkonium in a <C 1 cases. 
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Figure 3: Phase space portraits of the periodically driven: a) hydrogenlike atom ; b) quarkonium 
(a 3> 1 case); c) quarkonium (a C 1 case), for e cr /e = 0.5. 
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Figure 4: Phase space portraits of the periodically driven: a) hydrogenlike atom ; b) quarkonium 
(a 3> 1 case); c) quarkonium (a C 1 case), for e cr /e = 0.9. 
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